В.В. ЖДАНОВИЧ, В.И. ИБРАЕВ, С.Н. ЩЕКИН (ЗапСибНИИгеофизика)

К типу малоамплитудных поднятий относится Маслиховская структура на Утлорской площади (рис. 1). Целью работы являлось определение контура залежи пласта AC₁₁ структурно-литологического типа в коллекторах мелководно-шельфовой формации. Детальная съемка на участке в 96 км² проведена с применением методики регулярной площадной системы наблюдений.

На этапе интерпретации данных сейсморазведки главная задача – стратиграфическая привязка отраженных волн (ОВ), от правильного решения которой зависит точность последующих построений и достоверность геологических результатов. Стратификация ОВ возможна рядом способов, но наиболее точное решение дает система интерпретации комплекса геолого-геофизических данных, предназначенная для объяснения и прогнозирования волнового поля на основе создания эффективной геосейсмической модели реальной среды [2]. Система создания этой модели, реализованная в виде программных средств, входящих в комплекс прогнозирования геологического разреза [4], получила широкое распространение благодаря следующим факторам: 1) наличию расчетных значений плотностей (а) и пластовых скоростей (v_пл); 2) оценке влияния отдельных границ на формирование волнового поля, 3) установлению положения временного окна для погоризонтного динамического анализа.

Стратиграфическая привязка ОВ и определение физических свойств исследуемого объекта (пласт АС₁₁) проводились в двух точках по имеющимся данным ГИС в скв. 1 и 2, которые являлись единственными в период интерпретации. Благоприятным условием с позиции оценки свойств пласта АС₁₁ в этих скважинах по результатам модельных исследований является разнотипность насыщения. В скв. 1 пласт АС₁₁ представлен нефтенасыщенным коллектором, а в скв. 2 по данным ГИС – водонасыщенным, т. е. скв. 2 находится за контуром ВНК.

На рис. 1 показаны окончательные результаты моделирования для решения перечисленных выше задач. Совместное рассмотрение литолого-стратиграфической колонки и синтетических волновых полей, полученных разложением волнового поля на элементарные волны, а также исключением границ, дает однозначный результат стратиграфической привязки ОВ. В частности, исследуемый объект – пласт АС₁₁ в скв. 1 (границы 22 и 23), формирует отражение на времени экстремума 1,920 с с учетом фазовой поправки Dt_ф= 0,036 с. Сопоставление синтетического разреза с исключением границ и фрагмента временного разреза, а также их сходство подтверждает высокую степень достоверности подбора модели среды.

Результаты расчетов для уточнения петроакустических свойств пород на основе базовой модели приведены в таблице. Критерием оценки точности подбора базовых моделей варианта 1 в обоих случаях является коэффициент корреляции r₂, характеризующий степень подобия подобранной синтетической трассы и трассы, отображающей реальный геологический разрез. Данный единичный критерий используется при решении задач стратиграфической привязки, но статический подход оценки на основе ФВК не гарантирует достоверности расчетных значений петроакустических свойств пород. Поэтому оценочные расчеты проведены на основе многофакторного анализа с учетом времени экстремума отраженной волны (Т_э), коэффициента отражения на границе (К_отр), коэффициента корреляции между исходной синтетической и реальной трассами (r₁), коэффициента корреляции между подобранной синтетической и реальной трассами (r₂) и коэффициента пористости (K_п).

где s – плотность породы; d_т – минеральная плотность твердой фазы; d_ж– минеральная плотность жидкой фазы.

При изложенном выше методическом подходе с использованием программных средств пакета геосейсмического моделирования из комплекса СЦС-3 ПГР оценочным критерием достоверности является исходный коэффициент корреляции r₁. По колебаниям значений r₁ можно судить об изменении k_п.

При сопоставлении модельных значений параметров пласта АС₁₁ варианта 1 скв. 1 и 2 при равных прочих условиях отмечается уменьшение K_отр на кровле пласта в 2 раза и между значениями вариантов 5 и 2 скв. 1 и 2 в 3,6 раза. Пласт АС₁₁ в скв. 2 представлен водонасыщенным коллектором, а отражение от его кровли характеризуется пониженным значением К_отр при относительно высокой пористости пород. Совокупность данных факторов предполагает целесообразность прогноза зоны распространения продуктивного коллектора по распределению амплитудно-энергетических параметров с привлечением программных средств площадной динамической интерпретации, входящих в состав интерпретационной системы ИНТЕРСЕЙС [1].

Структурная карта по горизонту Д_в¹ (пласт AC₁₁) на рис. 2 характеризует приуроченность залежи к структурно-литологическому типу. Результаты ГИС и моделирования в разведочных скв. 1 и 2 дают основание полагать, что наибольший интерес для рассмотрения представляет контур пласта АС₁₁, ограниченный изогипсой –2190 м. О формировании и росте Маслиховской структуры в готерив-барремское время можно судить по структурной карте отражающего горизонта D_в¹. Так, на времени 1,95 с контур структуры уже вырисовывается, но в пределах изогипсы –2190 м динамически выраженные участки отсутствуют. Однозначно определяется область сноса осадков палеопотоками с юго-востока. На временах 1,94 и 1,93 с контур Маслиховской структуры определяется более отчетливо. Динамическая выраженность волнового поля в сводовой части Маслиховской структуры отмечается на времени 1,920 с. По результатам такого сопоставления возможно представить рост структур в разное геологическое время, но оценка полей динамических параметров остается при этом на качественном уровне.

Количественная оценка вариаций полей динамических характеристик проводится по картам, рассчитанным средствами ИНТЕРСЕЙС. Окончательным результатом совместной интерпретации совокупности двухмерных сейсмических полей в данном случае является карта площадного распределения комплексного параметра S_К (рис. 3), полученная как свертка одиночных величин (интервальная энергия, средневзвешенная частота, энергия в частотном диапазоне 8–32 Гц, скорость v_ПАК).

Распределение интервальной энергии по площади неравномерно, так в контуре, ограниченном изогипсой –2190 м, значения ее достигают уровня 30 усл. ед. Эта аномальная зона с небольшим пережимом по профилю 0–1841313 (пк 43, 75–48, 25) распространяется на запад, но уже на более низком гипсометрическом уровне. Частоты по отражающему горизонту Д_в¹ в пределах площади меняются от 36 до 48 Гц. Выделяемая низкочастотная аномалия на рассматриваемом участке площади занимает более обширную область в отличие от энергетических параметров и имеет значения 36–39 Гц. Распределение энергии в низкочастотной области идентично интервальной энергии со значениями от 3 до 33 усл. ед., а в рассматриваемом контуре аномальная зона ограничивается изолинией 12 усл. ед. Градиент скорости v_пак по результатам псевдоакустических преобразований (комплекс ПАК СЦС) составляет 100 м/с, аномалия v_пак в рассматриваемом контуре имеет пониженные значения от 2190 до 2950 м/с, что характеризует повышение истинного поглощения пласта АС₁₁ в сводовой части Маслиховской структуры. Карта распределения комплексного параметра S_К является результатом свертки рассмотренных выше одиночных динамических параметров. Вариация поля S_K представлена в виде центрированного распределения, а уровни меняются от –5 до 5 усл. ед. Результаты количественной оценки динамики волнового поля волны Д¹_В и данные геосейсмического моделирования дают основание полагать, что контур S_K, ограниченный изолинией 1, может быть идентифицирован с контуром нефтенасыщенного коллектора пласта АС₁₁, т. е. линией ВНК. Данное положение подтверждается бурением, проводимым по результатам выполненных исследований, являющихся составной частью работ по Маслиховскому месторождению. На рис. 4 представлена структурная карта по отражающему горизонту Д_в¹ (пласт АС₁₁) с результатами комплексной интерпретации и дополненная информацией по эксплуатационному бурению.

1. Интегрированная система комплексной интерпретации данных сейсморазведки и ГИС на ЭВМ / В.В. Жданович, А.П. Ковалев, Б.В. Монастырев и др. // Труды 33 Международного геофизического симпозиума.– Прага.– Т. 1.–1988.– С. 139– 147.
2. Метод эффективной сейсмической модели / Б.Я. Гельчинский, А.А. Белозеров, Н.И. Бердинникова и др.– Л.: Изд-во ЛГУ.–1975.
3. Петрофизическая характеристика осадочного покрова нефтегазоносных провинций СССР (справочник).– М.: Недра.– 1985.
4. Программно-методический комплекс для решения задач прогнозирования геологического разреза (СЦС-3-ПГР) / Г.Н. Гогоненков, А.Г. Авербух, С.Н. Птецов и др.– М.; Изд. ЦГЭ.–1984.

The solution of problems concerning stratigraphic confinement of reflected waves and the determination of physical propetries of the object under investigation is considered on the basis of creating an efficient geoseismic model. A quantitative evaluation is made for the fields of dynamic parameters in terms of defining the outline of pool in Neocomian terrigenous deposits at one of the Western Siberian fields. This paper contains an analysis of advantages to apply procedures from sets of programs "СЦС-3-ПГР" and INTERSEIS when processing and interpreting geologically the data on spatial seismic exploration.

1 – песчаники, 2 – глины, 3 – переслаивающиеся песчаники и глины, 4 – битуминозные аргиллиты, 5 – залежь нефти, 6 – кривая ПС, А – визуализация моделирования, В – разложение волнового поля на элементарные волны, С – синтетический разрез, последовательное исключение границ (2–51), D – фрагмент временного разреза

Рис. 2. Структурная карта по горизонту Д_в¹ (а), временной разрез (б) и горизонтальные срезы волнового поля (в):

1 – изогипсы по горизонту Д_в¹, м 2 – контуры наибольшей интенсивности динамических изохрон отражения Д_в¹

Рис. 3. Карта распределения комплексного параметра S_К горизонта Д_В¹ (пласт AС₁₁)

Рис. 4. Структурная карта Маслиховского месторождения по горизонту Д_в¹ (пласт АС₁₁)

1 – изогипсы отражающего горизонта Д_в (пласт АС₁₁), м скважины 2 – разведочная, 3 – эксплуатационная, 4 – аномально высокие амплитуды, выделенные при визуальном анализе сейсмической записи, 5 – повышенные значения комплексного параметра (динамический анализ), 6 – предполагаемый контур залежи пласта АС₁₁, 7 – принятый контур нефтеносности